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1. Introduction 

The European Twisted Mass Collaboration (ETMC) is generating gauge configurations with 
four dynamical flavors: two degenerate light quarks and a pair of strange and charm quarks with 
their masses approximately fixed to their physical values {Nf = 2+1 + 1). Several volumes and 
lattice spacings smaller than 0.1 fm are being produced [1, 2]. We present here first results using 
the A^j = 2 + 1 + 1 formulation to calculate observables probing nucleon structure. In particu- 
lar we show results for the non-singlet moment of the nucleon's unpolarized parton distribution 
{^)u-d- "^^^^ lowest non-trivial moment of the quark momentum distribution. There are no 
disconnected contributions, which are computationally very demanding, and hence can be 

computed relatively easily with lattice QCD. Furthermore, is also known accurately from 

global analyses of experimental measurements, thus it provides a good benchmark for lattice com- 
putations of nucleon structure. The successful determination of {x)u-d will provide confidence in 
the applicability of lattice QCD to predict other quantities of interest that might not be so well 
accessible experimentally. 

The operators necessary to calculate nucleon structure require a renormalization and matching 
to the continuum scheme used to analyse the experimental measurements. The renormalization is 
performed non-perturbatively in the RI'-MOM scheme and then matched perturbatively to the con- 
tinuum MS scheme to compare to the phenomenological results for the moments of parton distribu- 
tion functions. In particular, we discuss the non-perturbative determination of the renormalization 
factors of the vector and axial vector operators with zero and one derivative. 

Parton distribution functions (PDFs) are defined in Minkowski space and therefore are not 
directly accessible to lattice calculations. However, the moments of PDFs can be related to local 
operators that can be calculated in Euclidean space with lattice methods. The first moment of the 
parton distribution function, q{x,}JL^), is given by 

{^)q,li^= J dxxq{x,ix^) = J dx x{q{x,lJ.^) +q{x,lJ.^)] . (1.1) 

We consider the isovector combination (^)„_(/ in order to eliminate disconnected contributions. It 
can be determined by evaluating the expectation value of the local operator O^^ given by 

{p,s\qY^^iD''h'q\p,s) \,2 = 2{x) 2 p^^ p'l (1.2) 

Here q denotes the quark doublet {u,d) and is the Pauli matrix acting on the flavor indices. 

2. Lattice techniques 

We use the Wilson twisted mass fermion action, which has the advantage of leading to physical 
observables that are automatically ff{a) improved [3]. Our Nf = 2 + 1 + 1 setup, where we include 
dynamical up, down, strange and chaim quarks, maintains this automatic ff{a) improvement. Thus 
no operator improvement is necessary for O^^ and {x)^_j is accurate to ^(a^), thus providing an 
advantage compared to the improved Wilson actions that would require additional calculations for 
the operator improvement. In Table 1 we show the details of the ensembles used so far in this work. 
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The pion masses range from 320 to 450 MeV. The volumes satisfy nijiL > 4 in order to keep finite 
size effects small. The lattice spacing for the coupling j3 = 1 .95 has been estimated in the mesonic 
sector to be a ~ 0.078 fm. Further details are available in Ref. [1]. 
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Table 1: In the first column we give the bare coupling. The second column is the spatial extent of the lattice 
and the temporal length is always T = 2L. The bare quark mass in lattice units is given in the third column. 
The fourth column is the pion mass in lattice units with the statistical and systematic errors. The results of 
{x)u-d from this work are given in the fifth column and the sixth lists the number of configurations used in 
the calculation. 

For the lattice calculation of (x) we need the evaluation of the proton three-point correlation 
function given by 

1 + 74 



C^p (^sink ^source) ^operator) — Tr 



y~', {j p{tmn]^i^smk)0 (Xoperator)-^/) (-^source) ) ) (2-1) 

'^operator j^siiik 



where Jp (Jp) is an appropriate proton lattice creation (annihilation) operator. The conesponding 2- 
point function is given by an identical expression but with the omission of the operator insertion. A 
standard ratio of the three-point function to the two-point function gives the desired matrix element 



in the large Euclidean-time limit: ^operator — 4. 



Figure 1: C3p/C2p as a function of fopeiator 
for fixed source-sink separation fsink ~ fsource 
on the ensemble with mjt « 450 MeV. 
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The ratio of Csy, and C2p is computed as a func- 
tion of foperator and for large enough time separa- 
tions it becomes a constant yielding the matrix 
element of interest. This is illustrated in Fig. 1 , 
where we plot the dependence of the ratio on 
^operator " ^source for fixcd sourcc-sink Separation 
fsink — ^source = 12a. To evaluate Csp we construct 
sequential propagators at the sink. In this ap- 
proach one fixes the sink and source location, x 
and X4 respectively, as well as the quantum num- 
bers of the initial and final states. The sum over 
the spatial x' can then be done implicitly by per- 
forming an inversion with a source constructed 
using the forward propagators generated at the 
source and the state at the sink. 



The operator couples to the valence quai^ks at an intermediate time ^operator- Fixing the source- 
sink sepai^ation gives rise to contributions from excited states and, therefore, the source-sink sepa- 
ration has to be chosen sufficiently large to rule out those contributions but small enough to avoid 
the exponentially dropping signal-to-noise ratio. It is essential to use Gaussian smearing for the 
quarks in the proton interpolating fields Jp to increase the overlap of the creation and annihilation 
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operator with the ground state of the proton as well APE smearing for the gauge fields that enter 
the smearing functions to reduce noise. For the evaluation of the correlation functions we used the 
parallel contraction code ahmidas [4]. 



3. Results 

In order to assess cut-off effects we examine the 
dependence on the lattice spacing, a, of the nu- 
cleon mass which has been computed at three dif- 
ferent values of a. We show in Fig. 2 the nucleon 
mass as a function of (a/ro)-^, where ro is the 
Sommer parameter. More details can be found 
in Ref. [5]. The observation is that the i^(a^) 
scaling is mild and compatible with zero. Ad- 
ditionally the dependence of {x)^,_^ on the lattice 
spacing that we have found using the Nf = 2 en- 
sembles [6, 7] is negligible. Therefore, we have 
reason to believe that the lattice cut-off effects for 
{x)j^_j using the A^y = 2 + 1 + 1 will also be small. 
We are currently investigating this issue. 




Figure 2: The nucleon mass in units of ro 
as a function of (a/ro)^ for two reference 
pion masses. 



The renormalization constants are computed non-perturbatively in the RI'-MOM scheme at 
different renormalization scales using the momentum source method [8]. The advantage of this 
method is a high statistical accuracy and the evaluation of the vertex for any operator including 
extended operators at no significant additional computational cost. For the details of the non- 
perturbative renormalization see Ref. [9]. 

In the RI scheme the renormalization constants are defined in the chiral limit. Since the mass of 
the strange and charm quarks are fixed to their physical values in these simulations, extrapolation to 
the chiral limit is not possible. Therefore, in order to compute the renormalization constants needed 
to obtain physical observations, ETMC has generated Nf = 4 ensembles at similar lattice spacings 
so that the chiral limit can be taken [10]. We present the renormalization factors for the local vector 
and axial vector operators Zy and as well as for the one-derivative vector and axial vector 
operators, Z^y and Z^^, respectively. The later fall into different irreducible representations of the 
hypercubic group, depending on the choice of the external indices, jLt, v. Hence, we distinguish 
between Zdvi (Zdai) = -Z^y (^da) ^^^^ ^dv2 (Zdai) = Z^y^ (^da^)- Although we will use the 
Nf = 4 ensembles for the final determination of the renormalization constants it is interesting to 
compute the renormalization constants also in the A^/ = 2 + 1 + 1 theory and study their quark mass 
dependence. 

In Fig. 3 we show results for both the Nf = 4 and A^y = 2 + 1 + 1 ensembles in the RI'-MOM 
scheme. As can be seen, we obtain compatible values for all the operators. The same behavior is 
also observed in the case of Zy and Z^. This can be understood by examining the quark mass depen- 
dence of these results. In Fig. 4 we show, for the Nf = 4 case, the dependence of Zy, Za, Zdv, Zda 
on two light quark masses. The values we find ai^e consistent with each other. This explains the 
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Figure 3: One-derivative renormalization constants as a function of the momentum squared in lattice units 
in the RI'-MOM scheme. We denote Nf=4 results by the open circles and A^^ =2+1 + 1 by the crosses. 
Filled circles show the results after subtracting ff{a^) terms for the Nf = 4- case and x-symbols for the 
Nf = 2 + 1 + 1 case. 



fact that the results in the Nf =4 and N = 2 + \ + \ cases are compatible. Furthermore it makes 
any extrapolation of A^^ = 4 results to the chiral limit straight forward. 

The renormalization scale dependent Zdv and Zda need to be converted to the continuum 
MS-scheme, and for this we use a conversion factor computed in perturbation theory to &{g^). 
They are also evolved perturbatively to a reference scale, which is chosen to be (2 GeV)^. The 
constant that renormalizes {x)^^_^ is Zdvi- The results are shown in Fig. 5 both before subtracting 
the perturbative i^(a^)-terms and after. Using the subtracted results, the preliminary value of Zdvi 
in the MS-scheme at (2 GeV)^ is Zdvi = 0.998 calculated using the Nf = 4 ensemble with j3 = 1 .95 
and am^ = 0.194. 

In Fig. 6 we show the results for for several pion masses for the Nf = 2 + \ + \ en- 

sembles and give the values in Table 1. We compare with the Nf = 2 results for similar pion 
masses [11]. As can be seen, the Nf = 2 and A'y = 2 + 1 + 1 results are consistent showing that 
the effects of the strange and charm quarks in the sea are small compared to the, admittedly large, 
statistical eiTors. The results do not show a strong mass dependence and are higher than experi- 
ment. This behavior is consistent with the results of other collaborations [7, 12]. Thus calculations 
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Figure 4: Zy and Za (upper panel) and 
renormalization constants for the one- 
derivative operators (lower panel) as a 
function of the twisted mass. We have 
added 0.02 to the value of Zdai to dis- 
tinguish it from Zda2- 



Figure 5: One-derivative renormalization constants as a function 
of the momentum squared in lattice units converted to the MS- 
scheme for the Nf = 4 ensemble. Open circles denote results 
before perturbative subtraction of #"(0^) -terms and filled circles 
denote the i^(a^) -subtracted results. 



at even lower pion masses are needed to shed light on the behaviour closer to the physical point. 



4. Conclusions 

We presented results on {x)^_^ with Nf = 2 + 1 + 1 twisted mass fermions for pion masses in 
the range of 320 MeV to 450 MeV at j8 = 1.95 {a ~ 0.078 fm). We used non-perturbative renor- 
malization calculated for an A^y = 4 ensemble at the same coupling. The results are in agreement 
with those obtained with Nf = 2 twisted mass fermions. As expected from Nf = 2, there is no 
strong mass dependence and the results are higher than experiment for the pion masses used in 
these calculations. The chiral behaviour of {x)^_j will be studied using smaller pion masses and 
paying particular attention to excited state contributions in the three-point function. 
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Figure 6: {x)^_^ obtained from Nf — 2 and from Nj — 2 + \ + \ ensembles for a range of pion masses. The 
experimental value was taken from the ABKM09 analysis [13]. 
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